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ABSTRACT 



An unusually deep (V, /) imaging dataset for the Virgo supergiant M87 with the Hubble Space Telescope ACS successfully resolves its 
brightest red-giant stars, reaching M ; (//m) = -2.5. After assessing the photometric completeness and biasses, we use this material to 
estimate the metallicity distribution for the inner halo of M87, finding that the distribution is very broad and likely to peak near [m/H] 
=! -0.4 and perhaps higher. The shape of the MDF strongly resembles that of the inner halo for the nearby giant E galaxy NGC 5128. 
As a byproduct of our study, we also obtain a preliminary measurement of the distance to M87 with the TRGB (red-giant branch tip) 
method; the result is (m - M) = 31.12 ± 0.14 (d = 16.7 ± 0.9 Mpc). Averaging this result with three other recent techniques give a 
weighted mean rf(M87) = (16.4 + 0.5) Mpc. 

Key words. Galaxies: halos - Galaxies: individual - Galaxies: stellar content 



1. Introduction 

A key route to understanding the early enrichment history of 
large galaxies is through direct, star-by-star measurement of the 
metallicity distribution function (MDF) of its halo stellar popu- 
lation. Direct resolution of the old-halo stars with ground-based 
imaging is, however, generally restricted to the relatively small 
number of galaxies in the Local Group and not much beyond. 
The Hubble Space Telescope cameras with their superior spatial 
resolution and depth are able to extend this type of study consid- 
erably further, bringing many more galaxies of all types within 
reach. 

The rarest and biggest types of galaxies are the giant ellip- 
tical (gE) types, and to date, very few direct measurements of 
their MDFs have been made through their halo-star populations. 
Data of this type have been obtained with HST imaging only for 
the nearest gE galaxies, namely NGC 5128 at d = 3.8 + 0.1 Mpc 
(Harris et al. [lf99l 120001 [20021 120101 Rejku ba et a l. [2004b and 
NGC 3 379 atd= 10.5 ±0.8 Mpc (Gregg et al. 120051 Harris et al. 
2007b). To these we may add the MDF study of the intermediate- 
luminosity L eo ellip tical NGC 3377, at d - 10.8 ± 0.5 Mpc 
(Harris et al. I2007at . Mouhcine et al. d2005bl 120071) have used 
HST images to directly measure the halo-star MDF of nine spi- 
ral galaxies, finding a correlation between overall luminosity and 
the mean halo metallicity, and suggest that halos of massive spi- 
rals may not be dominated by metal weak stars. For this kind of 
work, the Virgo cluster of galaxies at d 16 Mpc is an especially 
tantalizing arena because it holds a rich collection of all types of 
galaxies including many large ellipticals. Among all of the Virgo 
members, however, M87 (NGC 4486) can be claimed to be the 
most intriguing target because it is the nearest example we have 
of a centrally dominant cD-type galaxy, which is likely accord- 



ing to the current galaxy formation framework (e.g. Thomas et 
al. 1200 5 ; Renzini 2006 for comprehensive discussions) to have 
undergone its major star-forming period very quickly in the early 
universe. 

In this paper, we take advantage of an unusual dataset in the 
HST Archive to probe the inner-halo stellar population of M87 
- that is, its fundamental population of old red-giant stars - for 
the first time. As will be seen in the discussion below, our abil- 
ity to gauge the MDF accurately is affected by the high level 
of photometric crowding in the images. Despite its limitations, 
however, it is good enough to permit some preliminary steps 
towards measuring the stellar MDF. In the following sections, 
we describe the raw data and our preprocessing; the photometric 
methods; measurement of the red-giant-branch tip and distance 
calibration; and finally a preliminary discussion of the metallic- 
ity distribution function. 



2. The Dataset and Preliminary Processing 

The raw images we use in this study are from HST program 
GO- 10543 (PI Baltz) and were taken with the ACS (Advanced 
Camera for Surveys) Wide Field Channel in the F6Q6W (wide 
V) and FM4W (wide I) filters. The 3.4' x 3.4' field, correspond- 
ing to 16 x 16 kpc linear scale at the Virgo distance, was cen- 
tered on M87. The exposures were taken in 61 separate space- 
craft visits, extending over a 71 -day period from 2005 Dec 24 
to 2006 March 5. These many exposures were originally de- 
signed to search for microlensing events. In F814W there are 
205 images totalling 73800 seconds, while in F606W there are 
49 images totalling 24500 seconds. Co-adding the individual ex- 
posures in each filter yields an unusually deep pair of images at 
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higher scale and a certain amount of gain in spatial resolution. 
Composite exposures of these images have already been used 
for three separate studies of the population of 2000 globular 
clusters that lie within the area of this ACS pointing (Waters et 
al. 120091 Peng et al. [20091 Madrid et al. [20091 . 

If the distance to M87 is 16 Mpc (see below), the expected 
magnitude of the TRGB (red-giant branch tip) at M/ = -4.05 
(Rizzi et al. 120071 1 is Itrcb — 27.0. As will be seen below, the 
limiting V,I magnitudes of these co-added exposures are eas- 
ily capable of directly resolving the TRGB and delineating the 
color-magnitude distribution (CMD) of the brightest RGB stars. 
The CMD can, in turn, be used to derive the MDF, as well as a 
direct calibration of the distance to M87 through the luminosity 
oftheRGBtip. 

For the purposes of this study, we carried out a new combi- 
nation of all the raw exposures in each filter through the APSIS 
software (Blakeslee et al. 120031 1. which performs accurate im- 
age registration, cosmic -ray rejection, and distortion correction 
with drizzle. Many hundreds of sources were used for measuring 
the offsets and small rotations of the individual exposures, yield- 
ing an accuracy of the registration of the raw images better than 
0.01 px. The unusually large numbers of exposures in each filter 
permitted subpixel resampling; we chose to construct final com- 
bined science images with a scale 07025 px 1 (half the native 
pixel size of the camera), and used the drizzle "point" (delta- 
function) interpolation kernel, which is equivalent to a square 
kernel with a pixfrac parameter of zero. Note that in the previ- 
ous globular cluster analyses of Peng et al. d20091 > and Madrid 
et al. ( 120091 1 we generated combined images with the Gaussian 
kernel and a scale of C'035 px . However, for the present study, 
we decided to rebuild the final science frames to gain the best 
possible effective resolution for faint-star photometry. 

A portion of our final science image in F814W, illustrating 
the degree of crowding in an outer corner of the field, is shown 
in Figure Q] The excised sample shown in the Figure is located 
in the southeast corner of the ACS field, at 12 30' 58.576" and 
+ 12°24'39.00", at a projected galactocentric distance of about 
150" or 12.1 kpc. The mean surface brightness of the region in 
FigureQ]in V and I bands is pv = 22.4 mag/sq arcsec and pj = 
21.2 mag/sq arcsec, from the surface photometry of Kormendy 
et al. 120091 (see also Zeilinger et al. 1 19931 >. The location of this 
field is shown in Figure [2] 

Another similarly deep probe of the RGB population in the 
Virgo region has been published by Williams et al. d20071 >. Their 
HST/ACS target field (from program GO- 101 31) is located in 
the intracluster region between M87 and M86, approximately 
180 kpc projected distance from each of these two giants. Their 
data consist of 26880 s exposure time in F814W and 63440 s 
in F606W, very similar total integration but with the opposite 
color emphasis than our central M87 field has. As Williams et 
al. discuss, the stellar population in this intracluster field is pre- 
dominantly metal-poor, and is likely to have a large mixture of 
material tidally stripped from many Virgo members. By contrast, 
the M87-centered field we analyze here should provide a rela- 
tively pure sample of the inner-halo RGB stars in this supergiant 
E alone and can be expected to be more metal-rich on average. 

3. Photometry 

To define the point-spread function (PSF) empirically for stel- 
lar photometry, we selected 22 moderately bright, isolated stars 
located across the ACS field and used these same stars for both 
filters (only about 40 bright stars in total are clearly distinguish- 
able in the entire field among the ~ 2000 globular clusters; see 



Madrid et al. 2009). Our target RGB halo stars are far fainter 
than these PSF stars, and the internal measurement uncertain- 
ties are dominated completely by sky noise and crowding, so 
the fine details of the PSF shape are not critical for this purpose. 
Nevertheless, we found no detectable differences in PSF size as a 
function of location in the field after the APSIS distortion correc- 
tions (see also Madrid et al. for further discussion). The FWHM 
of the PSF we obtained is 3.2 px (0.080") in F606W and 3.4 
px (0.085") in F814W, yielding a small but useful gain in effec- 
tive resolution compared with the nominal 0.1" resolution of the 
individual exposures. 

Across almost all the area of this M87-centered ACS field 
the level of stellar crowding is, to say the least, severe. It is only 
in the outermost corners that crowding is low enough to permit 
star-by-star identification and photometry by normal techniques. 
Once the PSF was established, we masked off all the area within 
3900 px (7.9 kpc or 97") projected galactocentric distance, es- 
sentially leaving only the four corners. This outer measured re- 
gion is marked in Fig. [2] To carry out the photometry, we used 
the IRAF daophot codes (see Stetson 1987) with a 2-pass se- 
quence of find/phot/allstar. For allstar the PSF fitting radius was 
adopted as 1 .5 px for the final runs to reduce crowding effects as 
much as possible, though we found that slightly larger values up 
to 2.5 px did not give significantly different results. Correction of 
the PSF-fitted instrumental magnitudes to large-aperture radius 
was done by curve-of-growth measurements from the same set 
of bright stars that were used to define the PSF. Lastly, these 
aperture-corrected instrumental magnitudes were converted to 
F6Q6W, F814W magnitudes and then into (V,I) with the cur- 
rent ACS/WFC filter zeropoints and conversion relations in the 
VEGAmag system (see Sirianni et al. 2003 and the ACS web 
pageffl), 

F606W = -2.51og(/ v ) + 26.421 (1) 
F8UW = -2.51og(/ 7 ) + 25.536 (2) 
V = F606W + 0.236(V - /) (3) 
/ = J F814W-0.002(V-/) (4) 

where fv,fi are the flux rates in counts per second for large- 
aperture photometry. The dominant uncertainty in the calibra- 
tion is in the correction from the small radii (2 px) that we are 
forced to use because of the high degree of crowding, to "infi- 
nite" large-aperture radius through the curves of growth. From 
these we believe the true external uncertainty in the photometric 
zeropoints to be +0.1 mag in each filter. 

Over the outer region defined above, our final data consist of 
105466 measured stars. The color-magnitude array for these is 
shown in Figure [3] subdivided into the four annular zones that 
are shown in Fig. [2] To facilitate comparison between zones, the 
expected level of the RGB tip luminosity for metal-poor ([Fe/H] 
< -0.7) RGB stars at Mj = -4.05 (Rizzi et al. l200"7l l is shown 
in each panel as the dashed line. For our present purposes, to 
minimize crowding effects as much as possible we restrict our 
following discussion to the two outermost zones (R gc = 115" - 
155"), containing 33890 measured stars. 

The RGB spans a broad color range, suggesting a large range 
in metallicity, or very large internal photometric scatter, or both. 
Our ability to see the true metallicity spread of the RGB tip is 
limited by the F606W exposures, which set the very distinct red- 
edge cutoff to the data at V < 30 (for similar cases, see the NGC 
3379 or NGC 5128 studies of Harris et al . 1200211200751 where the 

1 http://www.stsci.edu/hst/acs/analysis/zeropoints 
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Fig. 1. A 12.5" x 12.5" portion of the composite F814W image of the M87 halo, demonstrating its resolution into stars. Three 
globular clusters are seen in the lower half of the picture. This section is located in the extreme southeast corner of the ACS field, 
150" or 12.1 kpc from the center of M87. Stellar crowding increases rapidly toward the galaxy center but is manageable at this 
location. 




Fig. 2. Full ACS field of view. Lines corresponding to declination versus right ascension are drawn. Four annuli mark regions used 
in the photometry at 97" - 105", 105" - 1 15", 1 15" - 125", and 125" - 155", and are drawn around the center of M87. Note that 
the 155" circle does not show, because it is coincident with the bottom left corner of the ACS image. The central circle and the first 
two annuli are not used for final analysis due to crowding. Figure Q] is the portion marked by a square in the lower left corner. 



most metal-rich part of the population is cut off). A more specific 
demonstration of this point is shown in Figure[4] where the com- 
bined photometry for the two outermost annuli is plotted along 
with fiducial tracks for 12-Gyr-old RGB stars over the metallic - 
ity range [Fe/H] = -2.3 to +0.4. These tracks are the same ones 
us ed in previous studies of NGC 5 128 and NGC 3379 (Harris et 
al. 2002 2007b ) and are drawn primarily from the model library 
of VandenBerg et al. (120001 1 . At the right-hand edge of the CMD, 
the detection limits set by the V filter prevent us from measuring 
any stars more metal-rich than [m/H] -0.2 (i.e. stars in the 
range of the three reddest tracks in the model grid). In addition, 
at levels / > 27.6 the increasing photometric measurement un- 
certainties produce an asymmetric scatter in the observed colors: 
objects measured too blue by random errors fall well to the blue 
side of the most metal-poor model track; contrarily, any that are 
measured too red by similar amounts fall below the V limits and 



do not appear on the CMD. These arguments suggest that the 
MDF derivation needs to be restricted to only the very brightest 
section of the RGB to minimize biasses and excessive random 
spread. We will quantify this approach below after a brief dis- 
cussion of the TRGB and distance calibration. 

4. Determining the TRGB and Fiducial Distance 

Although the main focus of our M87 study is the metallicity dis- 
tribution, the data are also suitable for a preliminary estimate of 
the tip of the red-giant branch (TRGB). The bright-end termina- 
tion of the red giant branch has proven to be an extremely effec- 
tive standard candle for distance measurement to nearby galax- 
ies (well reviewed by Rizzi et al. 2007). Observationally it has 
the strong advantages that only single-epoch imaging is needed; 
that it is used for halo stars where internal differential redden- 
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Fig. 3. The color-magnitude distribution (7 vs. (V - /)) for the measured stars in the M87 inner halo. Data are plotted in four 
concentric annuli in projected galactocentric distance. These annuli cover a total radial range in R gc from 7.9 kpc to 12.5 kpc. The 
horizontal dashed line in each panel shows the predicted level of the RGB tip at M/ = -4.05 (see text). 



ing is usually not an issue; and that it rests soundly on a well 
defined and well understood breakpoint in the evolution of old, 
low-mass stars that are universally present in most galaxies. In 
addition, it is most easily employed in the / band where the level 
of the TRGB is nearly flat for the generous metallicity range 
[Fe/H] < -0.7 (see Rizzi et al. and papers cited there). 

The TRGB, which physically represents the helium-ignition 
point along the evolutionary tracks of the red-giant stars, is em- 
pirically defined by a sudden rise in the luminosity function 
(LF) of the halo stars. For recent examples of the sharp, near- 
ideal TRGBs that result from uncrowded photometry and sta- 
tistically large samples of stars, see McConnachie et al. (2005), 
Mouhcine et al. (2005a), or Harris et al. (12007b) . However, our 
M87 photometry is unavoidably in a different regime: because of 
the crowding and substantial star-by-star photometric uncertain- 
ties (which are already as large as +0.36 mag at the top of the 
RGB, as detailed in the next section), the TRGB does not appear 
here as a sharp onset in the LF but rather a smoothed-out rise. In 
Figure [5] we show the LF for the measured stars in the outer ra- 
dial range (1 15"- 155") of our M87 ACS field. Here, the data are 
smoothed with a Gaussian kernel of cr — 0.01 mag, though the 
result is insensitive to the particular smoothing width. The nu- 
merical first derivative (dn/dl) and second derivative (d 2 n/dl 2 ) 
or "edge response filter" ERF are shown below the smoothed LF. 

A prominent increase in the ERF amplitude first appears 
at / = 26.98 + 0.03 (internal uncertainty), which we identify 
as the TRGB. This apparent magnitude must be adjusted for a 
systematic photometric error (see next section) of Al = 0.12 
mag, yielding a true TRGB level IjTR GB) = 27.10. Adopting 
Mi(TRGB) = -4.05 + 0.10 (Rizzi et al. 120071 1. we obtain an ap- 



parent distance modulus (m—M)i = 31.15+0.14, where we have 
now included the +0.1 -mag uncertainty in the photometric ze- 
ropoint (see above). Subtracting Aj = 0.03 as listed in the NED 
then gives (m - M) = 31.12 + 0.14 and a true distance to M87 
ofd = (16.7 + 0.9) Mpc. 

As a secondary test of the degree to which crowding alone 
affects our result, we show in Figure|6]the smoothed LF and ERF 
for the four separate radial zones defined in Figs.|2]and[3] Aside 
from differences in the total numbers of stars, no significant off- 
sets are seen in the estimated TRGB from field to field. 

Given that the ERF technique is difficult to apply in such 
high crowding regimes, we estimated the TRGB level with an 
independent approach, similar to the ones used by Harris et al. 
(QMS! and McConnachie et al. d20"04l l. Here, the LF of the ob- 
served stars is modelled as the sum of an AGB and RGB for an 
adopted distance modulus, convolved with the photometric er- 
rors, and compared to the observed LF. The RGB is assumed to 
have a step function onset at the TRGB, rising exponentially to 
fainter magnitudes; the AGB is assumed to start one magnitude 
brighter than the TRGB, rising exponentially but with lower am- 
plitude than the RGB. (In practice, the AGB can be thought of 
as the sum of the true old- AGB stars plus blended pairs from the 
RGB.) The exponential slopes are set to A log 7Y/A7 = 0.96, 
and the normalisation of the AGB to the RGB is a free pa- 
rameter. We convolve this model by the photometric error and 
completeness as functions of / magnitude (determined from our 
tests with artificial stars), and compare it to the observed LF. 
We found good fits to the observed LF for an M87 TRGB at 
/ = 27.1 +0.2 and an AGB-to-RGB ratio at the TRGB of* 5-10 
percent. The AGB/RGB ratio plays only a minor role in the fLt- 
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(V I) 



Fig. 4. Color-magnitude distribution (7 vs. (V - I)) for the measured stars in the outer radial range (115" - 155") of our ACS field. 
The superimposed solid lines are fiducial 12-Gyr RGB tracks from the VandenBerg et al. 2000 model library, in steps of roughly 
0.1 dex from log (Z/Z ) = -2.0 to -0.1. These are supplemented with empirical tracks for two high-metallicity values at log (Z/Z ) 
= +0.16 and +0.4 (dashed lines; see Harris et al. 2002). 
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ting, which is most sensitive to I(TRGB). The fitting corrobo- 
rates I(TRGE) = 27.1 derived from the ERF method above. 

The stars lying above the RGB tip will be a combination 
of (a) field contamination from both foreground stars and faint, 
small background galaxies; (b) stars that are genuinely part of 
the M87 halo/bulge population, either young evolved stars or old 
AGB stars in highly evolved states including long-period vari- 
ables; and (c) accidental blends of two or more normal RGB 
stars. The field contamination component Nf is independent of 
the number Nrcb of M87 stars; the number of AGB stars goes 
in proportion to Nrgb\ and the number of blends varies as N^ CB . 
For our two outer zones, there are 761 stars in the 0.7-mag re- 
gion just above the TRGB, versus 4467 stars in the 0.5-mag re- 
gion just below it. For field contamination, comparison with the 



photometry of Williams et al. (120071) and Madrid et al. ( 120091 ) 
suggests that $ 40 foreground stars are present, a negligible 
contribution. The number of faint, small and misidentified back- 
ground galaxies is harder to assess without an adjacent control 
field to rely on, but comparison with the similar ACS ( V, I) Leo 
field data of Harris et al. (2007a) suggests very conservatively 
N ga i ^ 50, also negligible. 

For the NGC 3379 halo field, Harris et al. C2007bl ) found 
that the number of LPVs was 2.5% of the number of stars in the 
brightest 0.5 mag of the RGB. Scaling from that fraction then 
indicates we should expect Nlpv - 1 10 in our two outer zones. 
These plus the field contaminants then account for about a quar- 
ter of the CMD stars above the TRGB. There is no evidence for 
any young population of stars, which if present would produce 
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an upward extension of the blue side of the RGB (see Williams 
et al. i20071 for an example). 

Most or all of the remainder is likely to be due to blends: we 
expect (see Harris et al. l2007at 

Nl a 

N blend - (*) 

where N* is the number of stars in the field capable of gener- 
ating a blended pair brighter than the TRGB, a is the area of 
one resolution element, and A is the area of the field. We take 
A — 1 .42 arcmin 2 for the area of the outer two annuli enclosed 
within the field; a = 7.8 x 1(T 3 arcsec 2 ; and N* - 16000 for the 
number of stars between I = 27 and 28, the top magnitude of the 
RGB. These parameters predict about 200 blends lying above 
the TRGB. In fact this estimate is only a lower limit, because (as 
discussed in the next section) crowding-driven incompleteness 
of the photometry affects even the range / < 28 and N* could 
be quite a bit larger than our directly measured count. Though 
admittedly provisional, our conclusion is that the great majority 
of the objects brighter than the TRGB can be accounted for by 
the combination of the three normal factors given above. 

One useful comparison study of the distance to M87 is the 
photometry of the intracluster stars in Virgo by Williams et al. 
J2007I I. Although their target field is far from M87, it should in 
principle have a reasonably similar TRGB level given that M87 
is near the physical center of the Virgo environment. Although 
they do not specifically derive the TRGB, their CMD for the 
intracluster stars clearly shows I(TRGB) ^ 26.9 - 27.0, close to 
our value for M87 itself. 

Other indirect, but still quite relevant, TRGB-based calibra- 
tions of the Virgo distance include studies of a small number of 
its dwarf ellipticals. Harris et al. d 1998b obtained (m - m)o = 
30.98 + 0.20for VCC1 104 = IC3388, while Caldwell (2006]) ob- 
tained (m - M)o = 31.03 + 0.1 for seven other dwarfs. All of 
these are consistent with membership in the Virgo cluster, and 
with our result for M87. 

The only other direct distance calibrations for M87 that rely 
on resolved stars are (a) the observations of two novae (Pritchet 
& van den Bergh 1985), from which they obtain only the ex- 
tremely uncertain range of 30.4 - 32.4 for the distance modu- 
lus; and (b) the planetary nebula luminosity function (PNLF) by 
Ciardullo et al. ( 119981 1. After adjustment of their adopted zero- 
point M*(A5007) = -4.54 upward to M* = -4.67 to refl ect a 
change in the fiducial M31 distance (see Harris et al. 201^1, the 
PNLF modulus is (m - M)o = 30.92 + 0.16. The modulus is 
0.2 mag smaller than our TRGB value, though the two methods 
agree within their combined uncertainties. 

The largest and most precise compilation of relative dis- 
tances for the Virgo galaxies is by Blakeslee et al. (20091 see 
also Mei et al. 2007), from the surface brightness fluctuations 
(SBF) technique. SBF represents a "secondary" standard-candle 
method because it relies for its calibration on nearby galaxies 
whose distances are known beforehand from resolved-star meth- 
ods, although theoretical calibration via stellar population mod- 
els is also possible. Blakeslee et al. (2009) determine (m-M)o = 
31.11+ 0.08 for M87 in particular, starting from a mean distance 
modulus of 31.09 + 0.03 for the Virgo cluster as a whole. This 
in turn is based on the Cepheid calibration of the SBF method 
given by Tonry et al. (2001), corrected following Blakeslee et 
al. (|2002) to the final set of Key Project Cepheid distances 
(Freedman et al. l200TT) . 

Our TRGB distance measurement, and the earlier PNLF 
measurement, are gauges of the M87 distance that rely much 
less on Cepheids. The PNLF zeropoint, for example, relies more 



heavily on the fiducial distance to M31, for which close agree- 
ment now holds among several standard candles including RR 
Lyrae stars, TRGB, and Cepheids. Insofar as all these methods 
start from a common, fundamental Local Group distance scale 
that has been carefully assembled by a combination of many 
standard candles, the ~ +0.2-mag variances noted above are 
reflective of the small discrepancies that can emerge from the 
different methods in ways that are often still hard to pinpoint 
and that are not even consistent from one galaxy to another. For a 
more extensive discussion on another system (NGC 5 128) where 
several stellar candles can be accurately compared, see Harris et 
al. d20T0T >. 

Another direct, though somewhat less precise, comparison 
method of interest is the linear size distribution of globular clus- 
ters (GCs), which has been developed by Jordan et al. (2005) 
into a standard-ruler technique. The key quantity is the peak of 
the GC half-light radius distribution, normalized to host galaxy 
size and calibrated via the Milky Way GCs. Using the Jordan et 
al. data for M87 and their calibration, we obtain d = (16.4 + 2.3) 
Mpc or (m - M) = 31.07 + 0.30. 

Combining the four methods listed above (TRGB, PNLF, 
GC sizes, SBF/Cepheids), we obtain a weighted average dis- 
tance modulus (m - M) = 31.08 + 0.06 for M87, or d = 
(16.4 + 0.5) Mpc. A more precise TRGB distance especially 
could be obtained very straightforwardly with halo-star photom- 
etry in a less crowded region, and would in our view be the most 
effective way to calibrate the distance to this important galaxy. 

5. Tests of the Photometry 

To test the internal errors and completeness of the photometry 
we ran two separate artificial-star procedures. In the first series, 
stars were added to a representative 1500 x 500-px region of 
the image in the lower left corner (similar to the region shown in 
Figure 1). These added stars were distributed evenly in color and 
magnitude over the ranges 0.5 < (V — I) < 5.0 and 27.0 < I < 
29.6, as shown in Figure [7] These intervals deliberately covered 
a larger range in both color and magnitude than in our observed 
CMD (Figure [3}. The same daophot measurement sequence as 
done on the original frames was then carried out, with a two-pass 
sequence of find/phot/allstar. 

Figure [7] shows the results for inserted stars that were ac- 
tually recovered in the photometry, including both their input 
magnitudes and colors (center panel) and their actually mea- 
sured values (right panel). Of the total of 10000 added stars, 
just 1917 of these were successfully found and measured in 
both colors, and a high fraction of these lie in the upper left 
(bright, blue) part of the CMD. The completeness of detec- 
tion / = n{recovered) I niinput) as a function of magnitude is 
shown in Figure |8j the 50% completeness levels are reached at 
F606W = 29.00 and F814W = 28.15. To be classified as "re- 
covered" a star must be detected in both images. Note that at 
very faint levels (below the 50% point) the formal values of / 
tend to decrease rather slowly, a result of the very high degree of 
crowding. 

Figure (particularly the difference between the second and 
third panels) clearly indicates that internal random uncertain- 
ties are large at all levels of the CMD. Figure [9] displays the 
differences between the measured and input magnitudes more 
completely. For / > 28.0, one magnitude below the RGB tip, 
completeness of detection becomes low and the systematic er- 
rors increase, such that most stars are measured too bright. We 
do not consider this faint region further. The measured colors 
are on average slightly too blue independent of magnitude, by 
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Fig. 6. The luminosity function and edge response filter calculations for the four radial zones defined in Fig. [3]above. All have been 
determined in the same way as in the previous figure. 
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Fig. 7. Results of the first artificial-star set of tests. Left panel shows the color-magnitude distribution of the artificial stars inserted 
into the region of the image described in the text. Middle panel shows the stars from the left panel that were actually found and 
measured on both frames. Finally, the right panel shows the measured photometry for the stars that were recovered. The heavy 
dashed lines in the middle and right panels show the 50% detection completeness levels. 



(A(V - I)} = -0.22 mag. Over the top magnitude of the RGB 
(I - 21 - 28) the internal random scatter of the magnitudes and 
colors (rms) is <x = +0.36 mag in F814W, +0.53 mag in F6Q6W, 
and ±0.40 mag in (F606 W - 14 W). 

The artificially even distribution of input stars in Figure [7] 
overpopulates the brightest part of the luminosity function com- 
pared with the real CMD. A second run of artificial-star testing 
was therefore done, this time by inserting a set of RGB stars 
that followed a real CMD for the halo of a giant elliptical. We 
took the photometry for the halo stars of NGC 5128 measured 
by Rejkuba et al. (2005) with deep HST/ACS exposures that 
reach the horizontal branch and thus accurately populate the en- 
tire luminosity and color range of interest here. Of the more than 
77000 stars in the Rejkuba et al. dataset, about half (38640) were 
chosen randomly; these were shifted to the distance and red- 



dening of M87 to ensure that they occupied the same range in 
intrinsic luminosity and color; and were assigned random (x,y) 
locations in the same 1500 x 500 section of the M87 image as 
before. These added stars cloned from NGC 5128 increased the 
total light in the image section by less than 2%. 



The results of the remeasurement process are shown in 
Figure [10] Of 5200 input stars brighter than / = 29.8, just 570 
were recovered after the 2-pass photometry. Interestingly, the 
faintest ones recovered (7 > 29) often have very large system- 
atic errors that push them well upward in the CMD. However, 
the region / < 28 preserves much useful structure representing 
the metallicity range, losing only the reddest part of the range. 
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Fig. 8. The photometric completeness as determined from the artificial-star tests in the previous figure. The ratio / represents the 
fraction of input stars actually found and measured in both colors, and is plotted here as a function of magnitude in each filter. 
The curves drawn through the brighter parts of each set of points is a standard interpolation function of the form / = \[\ — 

a(m - mo)/ V 1 + a 2 (m - mo) 2 ] where mo is the 50% completeness level and a governs the steepness of the falloff. 
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Fig. 9. Random uncertainties in the photometry, determined from the first series of artificial-star tests described in the text. The first 
and second panels show the differences A/, AV (measured minus input) for the recovered stars, while the third panel shows the 
difference in color. 



6. The Metallicity Distribution 

The tests of the photometry described above lead us to restrict 
the MDF analysis to the brightest part of the RGB at 28.0 > / > 
27.0. In this range, even the random +0.40-mag scatter in (V—T) 
is much smaller than the 1.5-mag intrinsic spread in the actual 
CMD, enough to give us workable leverage on the underlying 
metallicity distribution. In addition, the =* 0.2-mag systematic 
error in the colors is far smaller than the intrinsic color range 
of the red-giant tracks as they fan out toward the top of the RGB 
(Figure|4]l, creating a shift at worst of 0. 1-0.2 dex in mean metal- 
licity. 

To estimate the combined effects of the quality of photome- 
try on the deduced MDF, we used the NGC 5128 data (Figure 



[Tol l again as a testbed. With the procedures described in de- 
tail in previous papers (Harris et al. 2002, 2007b, Rejkuba et 
al. 2005), interpolation within the grid of RGB tracks yields the 
MDF in the form of number of stars per unit interval in [m/H] 
= log(Z/Z ). The inner halo stars within 8 kpc in NGC 5128 are 
now well established to have a broad, predominantly metal-rich 
MDF peaking near [m/H] -0.3, and other giant ellipticals are 
expected to be similar (see the extensive discussions in Harris et 
al. 120021 Rejkuba et al. 2005 among others). 

The RGB-grid interpolation procedure applied to the NGC 
5128 stars is shown in three different versions in Figure [TJJ In 
the top panel, we take the 507 stars in Figure [Tol a) that lie in 
our adopted range 27 < / < 28 and derive the MDF in histogram 
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Fig. 10. Results of the second artificial-star set of tests. Left panel shows the color-magnitude distribution of the artificial stars 
inserted into the image section as described in Section 5. The stars are taken from the NGC 5128 halo field measured by Rejkuba et 
al. (|2005 ), shifted to the distance and reddening of M87. Middle panel shows the stars from the left panel that were actually found 
and measured on both frames; and the right panel shows the measured photometry for the stars that were recovered. 
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Fig. 11. Calibration tests of the metallicity distribution. In the top panel, the MDF is shown for the all the NGC 5128 stars from 
/ = 27 to 28 from the first panel of the previous figure. In the second panel, the MDF is shown for all the NGC 5128 stars that 
were recovered in the simulation tests (those in the middle panel of the previous figure). In the third panel, the MDF for these same 
recovered stars is shown, but now derived from their photometry as actually measured from the simulations (i.e. from the third panel 
in the previous figure). All three panels are normalized to the same number of stars. The bottom panel shows the ratio of the first 
and third histograms, normalized to an overall average of 1.00. 



form as shown there. This can be thought of as the intrinsic, input 
distribution free of photometric selection or scatter. 

In the second panel, the MDF derived through the interpo- 
lation grid is shown for only the 167 NGC 5128 stars that were 
recovered in the simulation tests (that is, those in Figure [TOl b)). 
The (/, V - I) values entered into the RGB grid are their true 
(input) values. The visible change in the MDF compared with 
the intrinsic version is the direct effect of the V-band cutoff, 



which is responsible for the loss of the highest-metallicity stars 
and thus also decreases the metallicity of the histogram peak. 

In the third panel, the MDF is now shown for the recovered 
stars at R gc = 1 15" - 155" or 9.3 - 12.5 kpc between / = 27 -28 
derived from their actually measured photometry (Figure fTOl c)). 
The overall histogram range and peak remain similar to the sec- 
ond panel, but relatively more stars appear at low metallicity 
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because of the systematic bias of the colors toward the blue as 
mentioned above. 

Lastly, the bottom panel of Figure QT| shows the ratio of the 
first and third histograms (input versus measured), normalized to 
an overall average of 1 .00. We use this ratio to make a reasonable 
correction of the MDF for both the incompleteness (which af- 
fects mainly the high-metallicity end) and the color bias (which 
affects mainly the low-metallicity end). 

The final results for our estimated MDF are shown in Figure 
[121 The 13972 M87 stars from Figure 0] in the range I = 27.0 - 
28.0 as described above are put through the interpolation grid 
of RGB tracks to construct the raw MDF shown in the bottom 
panel (shaded histogram). This distribution is then multiplied by 
the correction factor shown in the bottom panel of Figure [TT] to 
produce the fully corrected MDF shown as the heavy unshaded 
line. This corrected distribution represents our best estimate of 
the underlying MDF for the inner halo of M87, corrected for 
photometric incompleteness and color measurement bias. 

Even this corrected MDF can be relied on only for [m/H] 
< -0.3; any stars more metal-rich than that are almost totally 
cut off from our data and thus no correction factors are useful. 
At the opposite end, the MDF shape for [m/H] < -1.7 is very 
poorly determined and we can say only that there are very few 
stars present in that range. 

For direct comparison purposes we also show MDF data for 
the halo stars in three locations of the NGC 5128 halo, from 
Harris et al. ([19991 [20001 120021 . These were all taken with the 
HST WFPC2 camera in (V, /) and have rather similar limiting 
absolute magnitudes to our M87 ACS-based photometry, though 
they are less affected by crowding. The MDFs were derived in all 
cases with the same RGB grid of tracks and interpolation code. 
The mid-to-outer fields at projected galactocentric distances of 
21 and 31 kpc have MDFs that are virtually identical to each 
other and are combined in the upper panel of Figure Q/J] while 
the inner-halo 8 kpc field is shown by itself in the middle panel. 

Our corrected MDF for the M87 inner halo - at a mean pro- 
jected distance of 10 kpc - clearly resembles the 8-kpc NGC 
5128 field more closely than the outer fields. Over the range 
[m/H] < -0.3 where we can make the comparison, we conclude 
that the inner halos of both these giant ellipticals have basically 
similar MDFs that are broad, predominantly metal-rich, and with 
very small numbers of metal-poor stars. The MDF for M87 may 

reach a peak near [m/H] 0.4, but for the present this metal- 

licity should be viewed as a lower limit to the peak value. 

7. Conclusions 

We have used HST Archive data for an unusually deep set of 
F6Q6W, F814W ACS images to probe the red-giant stellar pop- 
ulation in the inner halo of M87. Although the crowding levels of 
the faint halo stars in this field are severe at best, the regions from 
115" - 155" (9.3 to 12.5 kpc projected galactocentric distance) 
give a useful first look at the brightest 1.5 magnitudes of the 
red-giant branch. This corresponds to distances of 1 AR e - 1 .9R e , 
where R e is the effective radius in the I band of 81" (6.3 kpc) 
(Zeilinger. [T993l . 

We have used this material to obtain a preliminary TRGB- 
calibrated distance to M87 for the first time. We find a distance 
d = (16.7 ± 0.9) Mpc, in good agreement with the few other rel- 
atively direct methods available, including the planetary nebula 
luminosity function, the Cepheid-calibrated SBF method, and 
the linear diameters of globular clusters. Photometry of a more 
outlying halo field than the one we are able to study here has 
the potential to give a more precise TRGB determination and 



would represent the most direct possible distance calibration to 
the Virgo members at the present time. 

As a second result of our study, we use the brightest 1 mag- 
nitude of the RGB stars to measure the metallicity distribution 
of the M87 inner halo for the first time. We find that it is very 
broad and predominantly metal-rich, with a peak that may lie 
very near the completeness limit of our V-band data but is at 
least as metal-rich as [m/H] -0.4. In general, the shape of the 
MDF strongly resembles the "8-kpc field" for the inner halo of 
NGC 5128, the nearby giant elliptical in the Centaurus group. In 
our view, these results strongly reinforce the view that the oldest 
stellar populations in giant E galaxies are fundamentally similar, 
regardless of whatever more recent accretion events may have 
happened to them. 

Our study and the Williams et al. (2007) photometry of an 
intracluster Virgo field, show that there is considerable poten- 
tial for extending this type of work to other Virgo members and 
the intracluster population as well. HST imaging with the ACS 
and WFC3 cameras is easily capable of resolving their halo red- 
giant populations and tracing the metallicity structures in this 
rich, nearby collection of large galaxies. 
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Fig. 12. Metallicity distribution for the halo stars in NGC 5 128 and M87. The upper two panels show the MDF (number of stars per 
0.1-dexmetallicity interval) from three locations in the halo of NGC 5128 (see text). In the bottom panel, the shaded histogram shows 
the MDF for the stars in the M87 inner halo in the magnitude range / = 27.0 - 28.0, uncorrected for photometric completeness or 
color bias. The heavy red solid line shows the MDF now fully corrected for those effects, while the heavy green dashed line shows 
the metallicity range in which the deduced MDF is most reliable. This corrected MDF for M87 is repeated in the middle panel 
(dashed line) for direct comparison with the inner-halo field in NGC 5128 (vertical scaling is arbitrary for display purposes). 
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